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Abstract 

The coupled dynamics of the scissors mode and the isovector giant quadrupole res- 
onance are studied using a generalized Wigner function moments method taking into 
account pair correlations. Equations of motion for angular momentum, quadrupole mo- 
ment and other relevant collective variables are derived on the basis of the time depen- 
dent Hartree-Fock-Bogoliubov equations. Analytical expressions for energy centroids and 
transitions probabilities are found for the harmonic oscillator model with the quadrupole- 
quadrupole residual interaction and monopole pairing force. Deformation dependences 
of energies and B(M1) values are correctly reproduced. The inclusion of pair correla- 
tions leads to a drastic improvement in the description of qualitative and quantitative 
characteristics of the scissors mode. 



PACS: 21.60.Ev, 21.60.Jz, 24.30. Cz 



1 Introduction 



An exhaustive analysis of the coupled dynamics of the scissors mode and the isovector giant 
quadrupole resonance in a model of harmonic oscillator with quadrupole-quadrupole residual 
interaction has been performed in pQ. The Wigner Function Moments (WFM) method was 
applied to derive the dynamical equations for angular momentum and quadrupole moment. 
Analytical expressions for energies, B(Ml)- and B(E2)- values, sum rules and flow patterns of 
both modes were found for arbitrary values of the deformation parameter. The subtle nature 
of the phenomenon and its peculiarities were clarified. 

Nevertheless, this description was not complete, because pairing was not taken into account. 
It is well known [2], that pairing is very important for the correct quantitative description of 
the scissors mode. Moreover, its role is crucial for an explanation of the empirically observed 
deformation dependence of E sc and B(M1) SC . 

The prediction of the scissors mode was inspired by the geometrical picture of a counterro- 
tating oscillation of the deformed proton density against the deformed neutron density [21 H] • 
Thus, as it is seen from its physical nature, the scissors mode can be observed only in deformed 
nuclei. Therefore, quite naturally, the question of the deformation dependence of its properties 
(for example, energy E sc and B(M1) SC value) arises. However, during the first years after its 
discovery in 156 Gd [5] "nearly all experimental data were limited to nuclei of about the same 
deformation (5 ~ 0.20 — 0.25), and the important aspect of orbital Ml strength dependence 
on 5 has not yet been examined", see ref. [B|. 

The first investigations of the 5-dependence of E sc and B(M1) SC were performed by W. 
Ziegler et al. [6], who have studied the chain of isotopes 148 > 150 > 152 > 154 Sm, and by H. H. Pitz et 
al. [7J and J. Margraf et al. [8], who have studied the chain of isotopes i 42 ' 146 - 148 ' 150 ]^. They 
found that the low-energy B(M1) strength exhibits approximately a quadratic dependence on 
the deformation S. 

Shortly afterwards it was discovered [9j [10], that in even-even nuclei the total low-energy 
magnetic dipole strength is closely related to the collective E2 strength of the 2f state and, 
thus, depends quadratically on the nuclear deformation parameter. 

Later J. Enders et al. [11] made a theoretical analysis of experimental data on the scissors 
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mode in nuclei with 140 < A < 200. Investigating the sum rules S + i and 5_i derived by E. 
Lipparini and S. Stringari [121 13], they found that the ratio Co = S + i/S-i is proportional to 
E sc with very good accuracy: E sc = 0.44a;. They also observed that the moment of inertia 
Jg S b of the ground state rotational band and the moment of inertia for the irrotational flow 
Juq = 5 2 J r ig (where J r i g is the rigid body moment if inertia) differ by nearly a constant factor 
(K ~ 10) over the entire region. Using this fact and identifying the giromagnetic ratio and the 
moment of inertia of the scissors mode with those of the ground state rotational band, they 
found with the help of the 5_i sum rule, that B(M1) SC is proportional to 5 2 . 

So, all the rather numerous experimental data demonstrate undoubtedly the S 2 dependence 
of B(M1) SC and the very weak deformation dependence of E sc . On the other hand, at the 
beginning of the nuclear scissors studies all theoretical models, starting from the first work by 
Suzuki and Rowe [14] , predicted a linear ^-dependence for both, B(M1) SC and E sc . 

It turned out that the correct 5-dependence is supplied by the pairing correlations. The 
effects of the pairing interaction in the description of the scissors mode were evaluated for the 
first time by Bes and Broglia [15J. They assumed "for simplicity that only the two subsets 
of levels which are closest to the Fermi level (n± and n± + 1) are affected by the pairing 
interactions". In this case B(M1) should be multiplied by the factor (u n± v n±+ i — u n±+ iv n± ) 2 ~ 
(u 2 l± — v 2 ± ) 2 = (e n± /E n± ) 2 with E { = J ef + A 2 and = — /i, where e« is a single particle 
energy, A is the gap and \i is the chemical potential. The value of (e n± /E n± ) was found by 
"making use of the fact that the moment of inertia is approximately 1/2 of the rigid body value 
obtained in the absence of pairing" . Thus, in accordance with the Inglis formula one has 

1/2 = J/J ng = [« _ v 2 J 2 /2E n± }/[l/2e n± } = (e n± /E n± f (1) 

and e n± / E n± =0.79. As a result B(M1) is reduced by the factor (e n± / E n± ) 2 = 0.62, i.e. the 
influence of pairing is quite remarkable. 

It was, however, noted by Hamamoto and Magnusson [16] that this result holds only for 
well deformed nuclei, where the equality ([T]) is valid. In general it is necessary to take into 
account the <5-dependence of the e n± / E n± - factor. This was done for the first time in ref . [16] . 
The authors applied "the method of averaging the position of the chemical potential between 
the occupied subshell (N,n±) and the empty shell (N,n± + 1)" to find that the ^-dependence 
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of B (Ml) is determined by the function 

Y = (5A A/3 ) (Vl + 4x 2 + — In IV1 + 4x 2 + 2x\) - 1 1 , 1 In \VTTx 2 + si | (2) 

V 2x / x \ xv 1 + a; 2 / 

with x = (Tu/JqS) I (2 A) . In the small deformation limit this function is proportional to S 2 , while 
for large 5 it deviates remarkably from such a simple dependence. The authors performed also 
a more realistic QRPA calculation for the Woods-Saxon potential with QQ and aa residual 
interactions, which confirmed their simplified analytical estimate. 

In [IT], N. Pietralla et al. established the 5-dependence of the e n± /E n± - factor phenomeno- 
logically. They were first to perform the theoretical analysis of the experimental data of the 
scissors mode in nuclei in the mass region 130 < A < 200. Following the idea of Bes and Broglia 
they parametrized the e n± /E n± - factor as 

E n Je n± = y/l + (bSy/{aS). 

The free parameters a and b were fixed by a fit to the experimental moments of inertia with 
the help of a formula equivalent to ([rj 

(finj_ / -^n± ) Jexpl ' Jrigi (3) 

where J exp = 3h 2 / E(2f) is the effective moment of inertia of the ground state band. In this 
way it was found that "the centers of gravity of the observed Ml strength distributions are 
always close to 3 MeV", i.e. "the data exhibit a weak dependence of the scissors mode on the 
deformation parameter" . They also derived a semiempirical formula for the total Ml strength 
of the scissors mode 

B(M1; 0+ - 1+) = 2.6c 2 rT |^-^ / 4 (4) 

(c g = 0.8 is the scaling factor of the giromagnetic ratio), which describes very well the exper- 
imental data and gives a deformation dependence "practically indistinguishable from the 5 2 
dependence" . 

A direct way to demonstrate the 5- dependence of the e n± /E n± - factor was suggested in 
[18J. E. Garrido et al. have shown that it is possible to extract analytically the 5 2 factor from 
the occupation coefficient $ a /3 = {u a vp — upv a ) 2 . Using the definitions of the u a , v a coefficients, 
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it is easy to write <& a p as 



e a - ep) 2 
AE a Ep 




(5) 



with 



P 



a/3 — 



(Eg — Epf 

(e a - ep) 2 



Uz-{z 2 -Ae 2 E 2 r 2 ] = e Ml + 




+ ...] 



(6) 



and E — \(e a — e ), e a p = ±(e a + ep),z = e 2 a/3 + E 2 + A 2 = \{e\ + e%) + A 2 = \{E 2 a + E 2 ). For 



the scissors mode e a — ep ~ food. The function (1 — P a p)/(E a Ep) has a regular dependence on 
S (no poles), so the coefficient $ a/ g and, respectively, B(M1) SC are obviously proportional to 5 2 . 
For the IVGQR e a — ep ~ 2^,o;, but its 5(M1) is proportional to 5 2 even without pairing due 
to other reasons (see section 3). 

In this paper we generalize the WFM method to take into account pair correlations. This 
allows us to obtain the correct (^-dependence for E sc and B(M1) SC in a slightly different way 
than in the papers, cited above. 

The paper is organized as follows. In section 2 the moments of Time Dependent Hartree- 
Fock-Bogoliubov (TDHFB) dynamical equations for normal and abnormal densities are cal- 
culated and adequate approximations are introduced to obtain the final set of six dynamical 
equations for the collective variables. In section 3 these equations are decoupled in the isoscalar 
and isovector sets and the isovector excitation energies and transitions probabilities are calcu- 
lated in the framework of the HO+QQ model. The results of calculations are discussed in 
section 4. Concluding remarks are contained in section 5. Some mathematical details are given 
in Appendices. 

2 Phase space moments of TDHFB equations 

The time dependent HFB equations in matrix formulation are [T9], [20] 



iMZ = [H, K) 



(7) 



with 
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The normal density matrix p and Hamiltonian h are hermitian; the abnormal density k and 
the pairing gap A are skew symmetric 

k ] = — k* , A f = -A*. 

The detailed form of the HFB equations is 

ihp = hp — ph — Ak 1 + kA\ —ihk = —hk — kh* + A — Ap* — pA, 

-ihp = h*p* - p*h* - A^k + tfA, -iTik = h*k ] + k ] h- A f + A^ + p*A f . (8) 

We will work with the Wigner transformation [20] of these equations. The relevant mathemati- 
cal details can be found in Appendix A. To make the formulae more transparent, in the following 
we will not specify the spin and isospin indices. The isospin indices will be re-introduced at 
the end. As a rule, we also will not write out the coordinate dependence (r, p) of all functions. 
The Wigner transform of jSJ) can be written as 

ihf = ih{h,f}- Ak* + kA*-— {A, + — {k,A*} 

2 2 

— [{{«, A*}} -{{A, «•}}] + ..., 

-ihf = ih{hJ}-A*K + K*A-^{A*,K} + j{K*,A} 

+y[{{K, A*}} -{{A, **}}] + ..., 
—itik = —hn — Kh — —{h,K,} — —{K,h} 

+a-a;-/a-|{/,a}-|{a,/} 

+y [{{h, *}} + {{«, h}} + {{A, /}} + {{/, A}}] + 
-ihk* = K*h + hn* + ^{K*,h} + l -^-{h,K*} 

-A* + /A* + A*/ + j {/, A*} + y{A*, /} 

-jKK. M) + {{h **}} + {{/, A*}} + {{A*, /}}] + (9) 

where the functions h, f, A, and k are the Wigner transforms of h, p, A, and k, respectively, 
/ (r, p) = /(r, — p), {/, g} is the Poisson bracket of the functions /(r, p) and g(r, p) (see 
Appendix A); the dots stand for terms proportional to higher powers of H. 
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To investigate collective modes described by these equations we apply the method of Wigner 
function moments. The idea of the method is based on the virial theorems by Chandrasekhar 
and Lebovitz [21] ; its detailed formulation can be found in J22J [23] . To study the quadrupole 
collective motion in axially symmetric nuclei it is necessary to calculate moments of Eqs. (JHJ) 
with the weight functions 



xz, p x p z , zp x + xp z , and zp x - xp z . 



(10) 



This procedure means that we refrain from seeking the whole density matrix and restrict our- 
selves to the knowledge of only several moments. Nevertheless this information turns out to 
be sufficient for a satisfactory description of various collective modes with quantum numbers 
K n = 1 + , as it was shown in our previous publications [TJ [22J E3]- In the case without pairing, 
this restricted information can be extracted from the TDHF equations and becomes exact only 
for the harmonic oscillator with multipole-multipole residual interactions. For more realistic 
models it becomes approximate even without pairing. The TDHFB equations Q are consider- 
ably more complicated than the TDHF ones, so additional approximations are necessary even 
for the simple model considered here. This is the subject of this section. 
Let us at first write out several useful relations: 

J d 3 p J d 3 rA{f, g} = -jd 3 pj d 3 rf{A, g} = - J d 3 p J d 3 rg{f, A}, 

Jd 3 pjd 3 rA{{f,g}} = Jd 3 pjd 3 rf{{A,g}} = J d 3 p J d 3 rg{{f, A}}, 

where A is any one of the above mentioned weight functions, / and g are arbitrary functions 
and {{/, g}} is defined in Appendix B. Integration of Eqs. (Q (including the terms of higher 
orders in ti) over the phase space with the weight A yields the following set of equations: 



ihj t Jd(p,r)Af = Jd(p,r) 



ih ~dl I f/(P ' V)A ' f / (/(p ' r) 



ih{A, h}f + A(A*k - k*A) - —({A, A*}k + {A, A}k* 

-y({{A,A*}}«-{{4A}K) 
ih{h, A}f + A(A*k — k*A) + —({A, A*}k + {A, A}k* 

-y({RA*}}K-{{A,A}K 




J d(p, r) A(h + 7i)k + ^{A, (h -h)} K - AA(1 -f-f) 




A}(/ - /) - \[{{A, (h + h)}}K + {{A, A}}(/ + /)] , 
* = J d(jp, r) -A(h + h) K * + j{A, (h - h)}K* + AA*(1 -f-f) 
+^{A,A*}(f-f) + ^[{{A, (h + &)}}«• + {{A, A*}}(/ + /)]1 , (11) 



where fd(p, r) = 2(2iTh)~ 3 fd 3 p Jd 3 r. It is necessary to note an essential point: there are no 
terms with higher powers of % in these equations. The infinite number of terms proportional 
to h n with n > 2 have disappeared after integration, as is demonstrated in Appendix B. This 
fact does not mean, that higher powers of h are not necessary for the exact solution of the 
problem. As it will be shown below, the set of equations ffTTl) contains terms, which couple 
with dynamical equations of higher order moments, which include, naturally, the higher powers 



It is convenient to rewrite the above equations in terms of h± = h ± h, f± = f ± /, 
A± = A ± A*, k ± = k ± K>*'. 



of h. 




+ 



J d{p, r) &+}/_ + {A, h-}f+) + A(A + k_ - K+ A_) 

4({{AA + }K-{UA_}K) , 



in— 
dt 




d{p, t)Ak + 




AA_(1-/+) 




ft 
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({{AMK + «AA-}}/ + ) 






ft 
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({{AMK + «4A + }}/ + ) 



(12) 
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These equations are strongly nonlinear, because A is a function of k (see, e.g., ref. [20J): 

d 3 p' 



A(r, p) = J (2 Jh)3 v (\ p ~ p 'l) K ( r ' p ')- 



(13) 



Having in mind small amplitude oscillations we will linearize: f± = f± + 8f±, k± = k° ± + 8k±, 
A± = A± + 8A±. The Hamiltonian should be divided into the ground state Hamiltonian h° 
and the residual interaction h 1 (and, if necessary, the external field). We consider h° without 
p-odd terms, hence h°_ = and as a consequence f°_ = 0. It is natural to take A real, i.e. 
A = 0, k°_ = 0. Linearizing (fT2l) and taking into account the last remarks we arrive at 



ift-n J d(p, r)A5f + = J d(p, r) 



ih 



({A, h°,}5f^ + {A, h l _}fl) + A(A°, 5k_ - k%8A. 



ih— J d(p, r) = J d(p, r 



■%.({{A,A° + }}5k--{{A,5A-}}kP + 
ih 



{{A,h° + }8f + + {Ah\}f + ) 



~({A,A° + }5k + + {A,5A + }kI) 



ih— J d(p,r)A5K + = Jd(p,r) 



ih 



Ah + 5n_ + y {A, hl}K° + - A5A_(1 - /}) 

h 2 



-{A, A° + }5/_ - —({{A, hl}}8K- + {{A, 5A_}}f + ] 



ih— J d(p, t)A8k- 
h 2 



d(p, r) Ah u + 5K+ + Ah\n° + - A5A + {1 - f + ) + AA u + 8f+ 



h» + }}8K + + {{A, fcUK + {{A, Al}}8f + + {{A, 8A + }}ti 



• (14) 



Until this point, our formulation is completely general. 

Now let us consider the popular case of pure monopole pairing. This means that the 
variation of the gap, 



5A(r,p) 



d 3 p' 



■v(|p-p'|)Mr,P') 



(15) 



(2vrft) 3 

will be projected on its monopole part. In the case of quadrupole vibrations, which we will 
study here, the variations 8f± and 8k± will have quadrupole multipolarities. As a consequence, 
when projecting formula (fT5"j) on the monopole part, the integral over angles will be equal to 
zero and we get 5A = 0. We also note that neglecting 8A corresponds to the usual Inglis 
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approximation. Then Eqs. (|14|) are reduced to 



ih 



{{A } h%}6f-. + {A,h 1 _}f?* 



n' 2 



thj t Jd(p,r)A5f_ 
d f 

i/i— J d(p, r)A5n^ 



ih— I d(p, r)A§n_ 



+AA° + 5k_ - —{{A,A° + }}5k 
d(p,r 



h° + w + + {a, hl}fl) - j{a, a° 



d(p,r) 



{A,h\}K° 



i l{A,A° + W---f{{A,h° + }}8K. 



d(p, r) 5k + + A/£k° + AA° 6f+ 

~({{A, h° + }}6K + + {{A, hl}}K° + + {{A, A° + }}5f, 



(16) 



To proceed further we are forced to do two approximations to get rid of higher rank moments 
and obtain a closed set of dynamical equations for second rank moments. First, the integrals 
/ d(p, r)Ah°,8K± contain fourth rank moments. The analysis of the integrand shows that we can 
neglect these integrals without a strong loss of accuracy. Indeed, the functions k± (and their 
variations) are sharply peaked at the Fermi surface, where the Hamiltonian h° + by definition is 
equal to zero. Therefore the product h? + 8k± should be small. Second, the (r, p)-dependence of 
A+ can generate, in principle, an infinite number of moments of various ranks. To simplify the 
problem we will consider here the commonly employed approximation of an (r, p)-independent 
gap A 5. = 2A = const, an approximation often used in nuclear physics and consistent with 
the monopole-monopole pairing force model. So, adding the isospin index r = (p,n), we finally 
have 



ihjJd(p,r)A6fl 

thjJd(p,r)A5r_ 
d r 

ih— / d(p, y)A8k t + 

(JjV J 



ih 



({A, h T + °}5fl + {A, /T 1 }/; ) + 2AA t 5k t _ 



ih— I d(p, r)A8n T _ 



d(p,r 
d(p,r 
d(p,r 

d(p,r) Ah?K? + 2AA T 5fl 



Til StO\ 



^({A,hl°}5fl + {A,h^}r + Q ) 
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~({{A,K?}}5kI + {{A,K*}}k?) 



(17) 



p 2 1 

We consider an axially symmetric model with h$ = h -m[c<j^ 2 (x 2 + y 2 ) + oj t2 z 2 ] — fx T 

and h\ = Z T (t)xz with Z n = XnnQ n + Xn P Q p and Z p = XppQ p + Xn P Q n , H T being the chemical 
potential of either protons (t=p) or neutrons (r=n). "Xtt 1 is the strength constant of the 
quadrupole-quadrupole residual interaction, Q T is a component of the quadrupole moment 
Q T (t) = Jd(p,r)xz5f T (r,p,t). It is supposed that \nn = Xw Obviously, in this model h T _ = 
and h T + = 2h T . Calculating the required Poisson brackets 

{xz, hi) = —L, {p x p z , h T } = -m(u T2 zp x + uj r2 xp z ), {zp x , h T Q } = —p x p z - mu T x 2 xz, 

{xp z , h T } = —p x p z - mu T z 2 xz, {xz, hi} = 0, {p x p z , hi} = -Z T (t)(zp z + xp x ), 
m 

{z Px , hi} = -Z T (t)z 2 , {x Pz , hi} = -Z T (t)x 2 , {{A, h T }} = 0, 

{{xz, hi}} = 0, {{ PxPz , hi}} = 2Z T (t), {{z Px , hi}} = 0, {{x Pz , hi}} = 0. 

we find, that the third equation of (ITT)) becomes trivial, giving four integrals of motion 
Jd(p,r)A5n+. Introducing the notation 

Q T (t) = Jd(p,r)xz5f T (r,p,t), Q T (t) = jd(p, r)xz5K T _(r, p, t), 

P T (t) = Jd(p,r)p x p z 5f T (r,p,t), P T (t) = Jd(p,r)p xPz 6K T Jr,p,t), 
L T (t) = J d(p,r)(zp x + xp z )5f T (r,p,t), L T {t) = J d(p, r){zp x + xp z )5K T _(r, p, t), 
Iy(t) = Jd(p,r)(zp x - xp z )5f T (r,p,t), I^(t) = J d(p, r)(zp x - xp z )5K T _(r, p, t) 
we find the following set of dynamical equations 



ihQ T 


= ih—L T + A T Q T , 
m 






~ T 

ihQ 


= AA T Q T + 2k T 4 Z T , 






%hP T 




l t - K 2 - 


-u: 2 )F y }+A^P\ 


ihp T 


h 2 

= AA T P T klZ T , 

2 






ihu 


= ih\— P T - miuol 2 - 


I - < 2 )Q T - 


((z 2 y + (x 2 y)z T 


ihi; 


= ihm{ul 2 - ul 2 )QT 


-th((z 2 y 


- {x 2 ) T )z T . 
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where 

(x 2 Y = J d(p, r)x 2 f°(r, p), (z 2 V = J d(p, r)^r°(r, p), 

fcj = J d{p, r)<(r, p), fcj = J d(p, r)xV<(r, p). (19) 

Eqs. ( 1181) will be simplified as far as possible to obtain results in analytical form. 

3 Simplified model 

The scissors mode is an isovector one, so it is natural to rewrite Eqs. (1181) in isoscalar and 
isovector terms. For the scissors mode, which we are interested in, the isovector set of equations 
can be decoupled from the isoscalar one with the help of the following approximations: 

A p ~ A n = A, k\ ~ kl = k A /2, ~ Aft = fco/2, 

(x 2 ) p ~ (x 2 ) n = (a; 2 )/2, (z 2 ) p ~ (z 2 ) n = (z 2 }/2, 5 P ~ 5 Q = 5, 

where 5 is the nucleus deformation. Introducing isovector variables Q = Q n — Q p , Q = Q n — Q p 
and so on, we can write the isovector set of equations as 

ihQ = ih—L + AQ, 
— m — 

ihQ = AAQ + 2k 4X iQ, 

ihP = -inm0 2 [(l + -)L-5Iy\ + AP, 

3 

h 2 

ihP = AAP--k oX iQ, 

2 2 & 

i%L = t h[-P-(2mQ 2 + -QooXi)(l + ^)Q}, 
m 6 6 — 

2 

ihly = -ih5(2mu 2 + -QooXi)Q, (20) 

o 

where Xi — \{Xnn ~ Xn P ) is the isovector strength constant. Usually one takes X i — a Xo> a 
being a fitting parameter. For the isoscalar strength constant Xo we will take the self consistent 
value 

3mtu 2 

Xo = —pr—. (21) 
voo 
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Following ref. [I] we take a = —2, i.e. a repulsive interaction with magnitude twice as large 
as the isoscalar one. Deriving (I20I) we used the self consistent expressions for the oscillator 
frequencies [1] 

u 2 x = u 2 y = u 2 {l + U\ ul = u\l- 2 -5\ co 2 = co 2 /(l+ 2 -S) (22) 

and the standard [22] definition of the deformation parameter 5 = 3<52o/(4<3oo), where Qoo = 
2 Jd(p, r)(x 2 + y 2 + z 2 )f°(r,p) = 2(x 2 ) + (z 2 ) and Q 2 o = 2(z 2 ) - 2(x 2 ) are monopole and 
quadrupole moments of nucleus, respectively. 

This set of equations has two integrals of motion : 

ih ~ rh 2 Xiko _ 2/1 W1 5., ifr 



2 m A £ + ll^T - M - + 3 "iAT^S " L = (23) 
and 

ly + 2mcj 2 5(l - a) — — — Q = consi. (24) 
— 4A + 2xi/c4 — 

Obviously these constants should be equal zero. By definition the variable Q is purely imaginary 
because k_ is the imaginary part of the pairing field k. Therefore Eq. (I24I) implies that the 
relative angular momentum I y oscillates in phase with the relative quadrupole moment Q of 
the imaginary part of the pairing field k. 

Analogously one can interpret Eq. ff23l) saying, that the variable L oscillates out of phase with 
the linear combination of two variables Q and P which describe the quadrupole deformation of 
the pairing field in coordinate and momentum spaces respectively. 

3.1 Eigenfrequencies 

Imposing the time evolution via e iQt for all variables one transforms Eqs. fTSOj) into a set of 
algebraic equations, whose determinant gives the eigenfrequencies of the system. We have 

Det = [E 2 - 2 X ihA){£ 2 - 2e 2 ) - 2e 2 (l - a)£ 2 - xi^h^/m 2 + 4ftV(l - a)5 2 = 0, (25) 

where S 2 = E 2 — 4A 2 , E = HQ, e 2 = h 2 u 2 {l + f). 

In the case of A = this equation is reduced to the known [Tj [25] equation for the scissors 
mode. In the case A ^ there are two solutions: 

E 2 ± = 4A 2 + [e 2 (2-a) + X ihA] 
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±\J[e 2 (2 - a) + Xih^} 2 ~ XiA[4fc 4 e 2 - k^/m?} - ATi 



4 cj 4 (1 - a)8 2 . 



(26) 



They describe the energies of the isovector GQR (E + ) and of the scissors mode (E_). 

It is worth noting that contrary to the case without pairing pQ the energy of the scissors mode 
does not go to zero for deformation 5 = 0. However this does not mean any contradiction with 
the known quantum mechanical statement that the rotation of spherical nuclei is impossible. 
It is easy to see from (1241) that the relative angular momentum I y is conserved in this case, 
I y = const, so the nature of this mode of a spherical nucleus has nothing in common with the 
vibration of angular momentum. The calculation of transition probabilities (see below) shows 
that this mode can be excited by an electric field and it is not excited by a magnetic field. Our 
estimation gives for the energy of this mode the value about 4 MeV, that agrees very well with 
the result of M. Matsuo et al [26J, who studied the isovector quadrupole response of 158 Sn in 
the framework of QRPA with Skyrme forces and found the proper resonance at ~4.2 MeV. 

It is known [TJ |2j that without pairing the scissors mode has a non zero value of an energy 
only due to the Fermi Surface Deformation (FSD). Let us investigate the role of FSD in the case 
with pairing. Omitting in (|20|) the variable P responsible for FSD and its dynamical equation 
we obtain the following characteristic equation 



Two solutions E 2 C = and Ef v = 4A 2 + 2«Xo^4^ + 2e 2 (l — a) reproduce the situation 
observed without pairing: the role of FSD is crucial for the scissors mode and is not very 
important for IVGQR. 

3.2 Transition probabilities 

The transition probabilities are calculated with the help of the linear response theory. The 
detailed description of its use within the framework of WFM method can be found in [lj, so 
we only present the final results. 

Electric quadrupole excitations are described by the operator 



E 2 [E 2 - 4A 2 - 2a X ohA - 2e 2 (l - a)] = 0. 



z 



(27) 



s=l 
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The transition probabilities are 

B(E2) U - 2| < HF 21 |0 > | - —— Qo0e ^ _ ^ 2 - £2(2 _ - - . (28) 

Magnetic dipole excitations are described by the operator 

z ~ eh 

For transition probabilities we have 

,~n, ,9 l-« 9^ ,9 E 2 .-4A 2 -2e 2 „ , . 

= 2| < ,K|0 > I = -ar-" _ ^ _ tH2 _ a) _ XikfA] A- (30) 

Multiplying B(M1) factors of both states by the proper energies and summing we find the 
following formula for the energy weighted sum rule 

— 9 

E SC B(M1) SC + E iv B{Ml) iv = (1 - a )-—Q 00 8 2 ^ N . (31) 

47T 

This expression coincides exactly with the respective sum rule calculated in [TJ without pairing. 
This means that there is no contribution to the sum rule which comes from pairing. This result 
can be explained by our approximation A = const. 

It is now a good place to discuss the deformation dependence of the energies and transition 
probabilities. First we recall the corresponding formulae without pairing: 




(E^r = 4rv 1 1 + 3 



8 ' M ' ! '^ tf ^) ^ (32) 

where the superscript "0" means the absence of pairing and we assumed a = —2. The scissors 
mode energy is proportional to 5, that becomes evident after expanding the square root: 



^ - M fi f 1 + H + mi + ■■■) • < 33 > 



where 5% = 1 + 5/3. 
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At a first superficial glance, the transition probability, as given by formula (1321) . has the 
desired (experimentally observed) quadratic deformation dependence. However, due to the 
linear ^-dependence of the factor E sc in the denominator, the resulting (^-dependence of B(M1)° SC 
turns out to be linear, too. The situation is changed radically when pairing is included. In this 
case the main contribution to the scissors mode energy comes from the pairing interaction (the 
term 4A 2 in fl26|) ). E sc is not proportional to 5 and the deformation dependence of B(M1) SC 
becomes quadratic in excellent agreement with QRPA calculations and experimental data [21 

SI Uni HBl 123 E3 EHl El] . 

The deformation dependence of B(Ml)i V is quadratic in 5, even without pairing, because 
the energy E iv is not proportional to 5 and depends only weakly on it. The pairing does not 
change this picture. 

4 Numerical results and discussion 

We have reproduced all experimentally observed qualitative features of the scissors mode. We 
understand that our model is too simplified to describe also precise quantitative experimental 
characteristics. Nevertheless we performed the calculations of energies and B(M1) factors to 
get at least an idea on the order of magnitude of the discrepancy with experimental data. The 
results of calculations for most nuclei, where this mode is observed, are presented in Table [I] 
and in Figures dHU Formulae (I26p and (1301) were used with the following values of parameters: 
a = -2, Qoo = A\R\ R = r Q A l '\ r = 1.2 fm, Q 2 = tu 2 /[(l + §5) 2 / 3 (l - fS) 1 / 3 ], hco = 41/ 'A 1 / 3 
MeV, h 2 /m = 41.803 MeV fm 2 . The gap A as well as the integrals fco, &4 were calculated by 
two methods: the semiclassical one, which is summarized in Appendix C, and the microscopical 
one described in Appendix D. 

Let us analyse at first the results obtained with semiclassical values of A, ko, k^ (columns 
th and / of Table 1). This method gives A « 1.32 MeV for all the nuclei considered here. 
The values of k , A; 4 vary smoothly from k = 37.7, k 4 = 4243.7 fm 4 for A=134 to k = 54.7, 
ki = 9946.3 fm 4 for A=196. Analysis of Table 1 shows, that overall agreement of theoretical 
results with experimental data is reasonable. It is, of course, not perfect, but the influence of 
pairing, especially on B(M1) values, is impressive. As it is seen, without pairing the calculated 
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Table 1: Scissors mode energies E sc (in MeV) and transition probailities B(M1) SC (in units of //fy); 
exp: experimental values, thl and th: full theory with ko, calculated by microscopical and semiclas- 
sical methods respectively, I: A ^ but ko = foj = 0, II: theory without pairing (i.e. A = ko = k^ = 0). 
The experimental values of E sc , 5 and B(M1) are from Ref. [T7] and references therein. El is explained 
in section 5. 
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energies (column II) are 1-2 MeV (1.5 - 2 times) smaller, than E exp , and B(M1) factors 
(column II) exceed experimental values 3-7 times. Taking into account A (without k&, ko) 
changes the results (columns I) drastically: the differences between calculated and experimental 
energies are reduced to 5 - 10% and calculated transitions probabilities are reduced by a factor 
of 1.5 - 2. Inclusion of (columns th) reduces the transition probabilities again by a factor of 
1.5, improving the agreement appreciably, and increases the energies by ~1 MeV deteriorating 
slightly the agreement with experimental data. The influence of ko is negligibly small, being 
of the order of ~ 1%. Finally we obtain that E t h exceeds E exp by ~1 - 1.3 MeV and that 
B(Ml)th exceeds B(Ml) exp approximately by a factor of 1.5 - 2. 

Let us analyse now the results obtained with microscopical values of A, k , fc 4 (columns thl 
of Table 1). Microscopical gaps for neutrons and protons vary in the range 0.8 - 1.0 MeV and 
1.0 - 1.2 MeV respectively, being appreciably smaller than the semiclassical ones. The values 
of &4 vary in the limits 4888 - 8936 fm 4 , being in the reasonable agreement with semiclassical 
values (in general they are smaller by ~ 10 - 20 %). Generally they increase with the atomic 
number A, however, due to shell effects it happens not monotonically, as in the semiclassical 
case, but quite irregularly. The integrals ko vary in the range 56 - 78. They are very sensitive 
to the detailes of level schemes - that is why they differ substantially (sometimes 60 - 70%) 
from the semiclassical values, being nevertheless of the same order of magnitude. 

The analysis of final results (Table 1) reveals the following trend: the scissors mode energies 
decrease (in comparison with those of column th) by ~ 0.5 - 0.7 MeV, that improves the 
agreement with experimental data, and B(M1) factors increase by ~ 10 - 20 %, that leads to 
a somewhat worse agreement with experimental data. 

What can be done to improve these results? The first step is obvious - it is necessary to get 
rid off approximations enumerated at the beginning of section 3, especially of the most crude 
one: A p = A n . As a result, it will be necessary to solve the coupled isoscalar and isovector 
sets of equations. The next possible step is to perform self consistent calculation with a more 
or less realistic interaction and taking into account the r-dependence of A. 

Another point, which should be clarified, is the role of the spin-orbit interaction. It is known 
[2], that experimentally observed low lying magnetic dipole strength consists of two separated 
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Figure 1: Scissors mode energies as a function of the mass number A for the nuclei listed in Table 1. 
The meaning of the symbols is explained in the caption of Table 1. 
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Figure 2: Scissors mode energies as a function of the deformation 5 for the nuclei listed in Table 1. 
The meaning of the symbols is explained in the caption of Table 1. 
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Figure 3: Scissors mode transition probabilities B(Ml) as a function of the mass number A for the 
nuclei listed in Table 1. The meaning of the symbols is explained in the caption of Table 1. 
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Figure 4: Scissors mode transition probabilities, normalized by a factor of A 4 / 3 , as a function of the 
deformation 5 for the nuclei listed in Table 1. The meaning of the symbols is explained in the caption 
of Table 1. 



20 



parts: orbital excitations in an energy interval ~2 - 4 MeV and the spin-flip resonance ranging 
from 5 to 10 MeV excitation energy. So, for the full description of the scissors mode dynamics 
it would be necessary to consider also the spin degrees of freedom. One may expect that 
the orbital part of the Ml strength (scissors mode) will be pushed down by the spin-orbit 
interaction, in agreement with experimental data. 

5 Currents 

According pQ the components of infinitesimal displacements in the plane y = are given by 

Ut) = V2BQ(t)z, Ut) = V2AQ(t)x (34) 

with 

A=^[l-2^(l-a)5]/[goo(l-|*)], 

B = A=[1 + 2^(1 - a)5]/[Q 00 (l + p)]. (35) 

It is useful to write these displacement as the superposition of a rotational component with the 
coefficient a and an irrotational one with the coefficient b 

— * 

£ = ae y x r + bV(xz) = a(z, 0, —x) + b(z, 0, x) — ► £ x = [b + a)z, £ z = (b — a)x. (36) 

Comparison of fl34l) with fl36l) gives 

b + a = V2BQ, b-a = V2AQ — ► b = (B + A)Q/y/2, a = (B - A)Qjy/2. 

Using here expressions (1351) one finds 

6 = 7(^3-^), a = 7 5(l-6 3 /g), (37) 

where g = E 2 /(6h 2 u 2 ), 5 3 = 1 + 5/3 and 7 = 3Q/[O 00 (l - + §5)]. The coefficients a, b are 
shown schematically as the functions of g on the figure 5. In spite of the simple behaviour of two 
curves, they describe rather interesting phenomena. There are two critical points: g\ = 5 2 /$3, 
where 6 = 0, and g% =83, where a = 0. They divide g axes into three regions. 
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Figure 5: Amplitudes a and b (in arbitrary units) of the rotational and irrotational current components 
as the functions of energy (g = E 2 /(6h 2 ui 2 )). 
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In the region < g < gx the rotational component of the current dominates. The ratio b/a 
is changed from 6/63 at g — to zero at g = gx- At this point the motion is pure rotational! 
The energy corresponding to this point is E\ = TuJjSJft/ 8%. Its value for various nuclei is shown 
in the Table 1. It is interesting that all experimental values of E sc are disposed just in the 
vicinity of Ex, the energies of nuclei with the small deformation (5 < .19) being disposed on 
the right hand side of this point, while the energies of nuclei with large deformations (5 > .2) 
are disposed on the left hand side. The theoretical value of E sc obtained without pairing (see 
formula ( 1331) ) is E® c ~ \Ex, i.e. it is disposed far to the left of Ex for all nuclei independently 
of their deformations. Inclusion of pairing shifts the scissors mode energy E sc strongly to the 
right so, that it becomes larger than the critical value Ex, but quite close to it, especially in the 
well deformed nuclei. One can note, that the fact of the proximity of theoretical and, especially, 
experimental values of E sc to the critical point Ex gives some grounds for the two rotors model 
of the scissors mode [I], where the pure rotational motion of the nuclear matter is supposed. 

In the region gx < g < g 2 the rotational and irrotational components of motion compete 
on an equal footing - the dominance of the rotational motion at the point g — gx is gradually 
replaced by the dominance of the irrotational motion at the point g = g 2 , the strengths of both 
components being equal ( i.e. \a\ = \b\ ) at the point g = 5. The energy corresponding to 
the point g 2 is E 2 = y/65~3?iLU = ^/^5l\/2huj. The energy centroid of IVGQR (without pairing) 
is E iv = 2\/2huj. For normally deformed nuclei (5 < 1) the factor ^35^ is smaller than 2. 
Therefore E iv is disposed on the right hand side of the point g 2 , i. e. in the region g 2 < g < 00, 
where the irrotational component of the current dominates. The inclusion of pairing shifts E iv 
more to the right. With g increasing the coefficients a and b grow gradually aspiring to their 
asymptotical values j8 and 7^3 respectively. 

One more detail of the interrelation of the two components of the current. Let us observe 
the motion of the ends (tips) of the scissors. Their displacements are determined essentially 
by the value of £ x = (b + a)z. In the region < g < gx coefficients a and b have the opposite 
signs. This means that two components of current move out of phase - the shear (irrotational) 
component bends the bodies of the scissors trying to resist to their rotation. This process is 
described very well and illustrated by the figures 1, 2 in the paper by R. Hilton [3]. In the 
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region gi < g < g 2 coefficients a and b have the same sign. Therefore, both components move in 
phase - now the shear component bends the bodies of the scissors reinforcing their rotation. In 
the region g 2 < g < oo coefficients a and b again have opposite signs and the situation returns 
to the case of the first region. However, now the shear component dominates, hence it is better 
to say that the rotation trys to resist to the shear. 

6 Concluding remarks 

The low energy magnetic dipole strength produced in all QRPA calculations (even in the 
schematic model HO+QQ) is always distributed over several states in the region < E < 4 
MeV, because each af3 pair contributing to B(M1) SC occurs at a different energy E a p = E a + Ep 
(see Introduction). The same picture is observed experimentally. These facts were in sharp 
contradiction with the anticipated single peak, predicted in early papers [HI [TJl US] and caused 
"controversy about the collectivity and the correct interpretation - in the sense of classical 
motion - of the orbital magnetic dipole modes [28]". It became clear that one needs some kind 
of a bridge between the classical and quantum mechanical approaches. For the description of 
this phenomenon "a formalism is required in which the interplay between the collective and 
single particle aspects of the system are adequately treated [29]" . 

Some authors "have tried to complement RPA calculations with realistic forces by classical 
and semiclassical methods which are flexible enough also to admit other solutions than the 
one anticipated [28]". This way however does not give a satisfactory solution of the problem, 
because usually the phenomenological models have not the direct connection with quantum 
mechanical ones. The most popular way to extract the underlying physical nature of the state 
from microscopic (RPA or QRPA) calculations is the calculation of overlaps [2J. However, the 
overlaps usually do not give the full information about the various properties of the studied 
phenomenon. As a matter of fact, this procedure answers only the question: "What part of the 
mode strength is excited by this operator (with which the overlap is calculated)?" Usually one 
calculates the overlap with the "synthetic scissors" [H [25] . It is known, however, that due to 
the coupling with the IVGQR [12] , the low-energy magnetic dipole mode contains rather big 
isovector plus isoscalar admixtures of two orthogonal shears, as shown in [29J by calculating 
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the overlap with the proper operator. Nevertheless, this overlap does not exhaust all strength. 
What other operators are necessary to clarify completely the nature of the considered mode? 

We think that the combination of two complementary methods, namely RPA (or QRPA) 
and WFM method, can be very useful. Based on the same approach (time-dependent HF 
or HFB together with a small amplitude approximation), they allow one, starting from the 
same Hamiltonian with the same forces, to obtain either (by RPA calculations) the refined 
microscopic structure or (applying WFM method) the crude "macroscopic" picture (which 
reminds very much the results of semiclassical approaches) of the same phenomenon. This 
interrelation of RPA and WFM was investigated in [30]. In particular, the identity of both 
methods in the case of a schematic model was demonstrated there. With the help of the WFM 
method, taking into account moments of higher and higher rank, one can produce a more and 
more detailed description of the phenomenon, achieving (at least in principle) the maximally 
fragmented picture given by the experiment (and RPA). 

One more remark is in order. Discussing the scissors mode energy and its ^-dependence, 
one has usually in mind the mean excitation energy (centroid) defined as the center of gravity 
of the Ml strength distributed among the low-lying K 71 = 1 + states 

E sc = ]T EiB(Ml)i/ B(Ml)i. 

i i 

The sums are evaluated in the energy intervals around 3 MeV. The evaluation interval is not 
strictly determined. For example, Hamamoto and Nazarewicz [25] took < E < 10 MeV 
in their calculations of superdeformed nuclei. Applying the WFM method, which produces 
centroids of resonances, we avoid such problems. 

In conclusion, the WFM method is generalized to take into account pair correlations and 
is used to calculate energies and transitions probabilities of the scissors mode. Excellent quali- 
tative and reasonable quantitative agreement with experimental data is obtained. In addition 
the interrelation of microscopic and semiclassical features of the scissors mode is clarified. 
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Appendix A 

We remind that each ingredient of Eqs. (jHJ) is a matrix in coordinates and spin projections. For 
example, k = k Tj(r y )(7 i. The spin structures of the normal and abnormal density are explicitly 

(Pt,t> \ / k r y\ 

V p r y J \-Kr,r> / 

or 

Pr,T;r',T = Pr,|;r',|, ^r,T;r',| = ~ «r,J,;r',T " 

With the help of the standard expression for the product of two matrices 

{AB) Tj(r)T i )ff i = J d 3 r" ^2 4,qr",o"5r'y' ; ry = / d 3 r"[A r ^yi ^B r » ^y ^> + A Ta . r r>^B r r/^.y (T t] 

cr" 

we find 

(^p)r,T;r',T = J ^ 3r "^r,T;r",TPr",T;r',T 5 (Mr.tirU = (A/$) r 4 ;r / )T = 0, 
(hp)r,l;r',l = / h r ,l;r",lPr",l;r',l, 



(Ak) r ^y^ = J d 3 r" A r ^. r >/^k r /i^.y^, (Ak) r ^y^ = (Ak) r ^y^ = 

(^)r,|;r',| = / dV' 'A P) j,. r // ^k r ii fry 4. 

As an example we write out the first equation of (jSJ) for a =] in detail: 



^Pr,T;r',T j d T /i r ,T;r",T/V',T;r',t — Pr,T;r",T^r",T;r',T 

A r ,T;r",J.« r //4 ;r / ) -r + £r,T;r",|Aj!„ ii;r , jT . (38) 

The Wigner Transform (WT) of the single particle operator matrix F Tlj(T - iT2 y is defined as 

[•^r 1 ,o-;r 3 ,cr / ]wT = ^ct.ct' ( r ; P) = J d 3 Se lP S ^ F r+s / 2 ,a;r-s/2,a' 

with r = (ri + i"2)/2 an d s — r i — r 2- It is easy to derive a pair of useful relations. The first 
one is 

K,A r >P) = / d 3 se^P: +s/2 ^_ s/2 ^ = J d 3 se-^ h F:„ s/2 ^ r+s/2 ^ 
= J d 3 se~ ip ' s/h F r t +s/2 a ,. r _ s/2 a = [F^ 1>(T/ . r2iCT \ WT 
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Le - [Fl^a;r 2 ,o>]wT = [F Tl ,*>;r 2> a\ * W t = F*, a {r,p). The second relation is 

iv CT ,(r,p) = F^(r,-p) = J d 3 se ip - s/h F r+s/2 ^ s/2 , a , 



r+s/2,cr';r-s/2,crJ ' 



For the hermitian operators p and h this latter relation gives [p* a )wT = Paa( r , — p) and 
The Wigner transform of the product of two matrices F and G is 



i 



h -A 



[FGVt = F(r, p) exp ^- A j G(v, p), (39) 



where the symbol A stands for the Poisson bracket operator 



^ ^dr-j dpi dpi dri 
For example the Wigner transform of Eq. (I38I) up to linear order in h is 

ift/r,r( r >P) = ^{ /l T,T( r .P)>/T,T( r 'P)} 

-A Ta (r, p)«| >4 (r, p) - — {A u (r, p), «* 4 (r, p)} 

+K T)i (r, p)A^(r, p) + y{K U (r, p), A* a (r, p)}, (40) 

where {/,<?} = / A g = f tt~tt' — Tr~Tri is the Poisson bracket of arbitrary functions 

^Xdridpi opidrij 

/(r, p) and <?(r, p); /i(r, p), /(r, p), A(r,p) and ft(r, p) are Wigner transforms of h ri , r2 , Pn,r 2 , 
A ri r2 and «; ri)r2 respectively. The functions ft, and / are real, because the matrices h and p are 
hermitian. This example demonstrates in an obvious way that the dynamical equations (JSj) for 
the matrix elements p r ,o-;r',<7; Ptvr' <?■> K r,±a;r',^a, K £ ±<rr' =p<t> w ^ n a = T an d °~ = l are transformed 
into eight dynamical equations for their Wigner transforms: 4 equations for / CT)CT (r, p), fa,a{*, p), 
/^-^(r, p), K*_ CT (r,p) with cr =f and 4 equations with a =|. By definition /^(r, p) = 
f a ,a( r , — p)- In the absence of spin dependent forces both of these subsets coincide and we can 
consider any one of them. 
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Appendix B 



Let us consider in detail the arbitrary term of equations ( TlTl) proportional to h n . Such terms 
appear after expanding the exponent in formula ( 1391) . After integration with an arbitrary 
function A(r, p) we have 



In 



d(p,r)A[fA g}= ld(p,r)[A 
dA 



df 9g 

dxi dpi 



df dg 

dpi dxi 



d(p,r) 



dpi 



Of ~n-l 
OX,; 



<9A 



9/ -"-I 
OPi 



A 

= ^4p + As- 



(41) 



Repeating the integration by parts we find 

<9 2 / ^«-2 <9g 



d(p,r 
d(p,r) f 



OA 
dpi 
d 2 A 
dpidpj 



dxidxj 



A 



<9p 



2 f <-»n-2 Qn 
A 



<->n— 2 

A s 



-4, 



d(p,r) 



dxidpj 



d 2 f 

dxidxj 

d 2 f dg 
— A — 

dpidxj dpj 

d 2 f 

' A 



dxidpj 



d 2 A 



dx^ 



dpidxj 
d 2 f 



d 2 f 



*n-2 



dxidpj 

<->ra-2 dg 



A 



+ A. 



dpidx 



dpidpj 
d 2 A 



A 



dxidxj 



dxj j 
A g 



dpidpj 



A xp + A xx . (42) 



Repeating once again the integration by parts we get 



a 3 / 



dxSxSx 



A 



-3 £>0 



- / d(p,r) 



<9 3 A 



dpidpjdpk 



d 3 f 



iUXjUXk 

A g 



dpk_ 
d 3 A 



d*f 



dxidxjdpu 
d 3 f 



A 



-3 90 



dpidpjdxk 



dxidxjdpu 



dx k 

«-^n— 3 

A 



(43) 



dxidxjdx^ 

It is easy to see from the structure of A pp (and A px , A xp , A XX) ) that in the case, when A(r, p) is 
a polynomial of an order k, all integrals I n with n > k are equal to zero. In our case k = 2 and 
we find from the above formulae, that ii = fd(p,r)A{f,g}, I2 = fd(p,r)A{{f,g}}, where 
{/, 0} is defined in Appendix A and 

~2 3 / Q2f Q2 Q2f Q2 Q2 f Q2 \ 



i ~ 1 \dridrj dpidpj dridpj dpidrj dpidpj dridrj 



Appendix C: Thomas- Fermi approach to nuclear pairing 



Following Ref.[3T], we define the density matrix averaged on the energy shell as 
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h = jr^-ME -i?)~E^- (H- (44) 
g(E) g\E) v 

which is a smooth function of E since 5 denotes a smeared delta function. The smeared level 

density g(E) (per spin and isospin in this paper) in the denominator of expression (|44j) ensures 

the right normalization of pE- The smooth quantities entering in 044[) are evaluated by replacing 

H, the independent-particle Hamiltonian, by its classical counterpart H c \ which corresponds to 

the Thomas- Fermi (TF) approximation [201 [32]. This approach is not limited to the evaluation 

of expectation values of single particle operators. Also the average behavior of two-body matrix 

elements can be calculated [31]. In this paper we are interested in the semiclassical evaluation 

of the average pairing matrix elements which at TF level read 

v(E, E') = ~ {E y g{EI) E HE - e v ) ~6(E' - e u ,) v)\v\${v\ u')), (45) 

where |$(i/, z/)) is an antisymmetric normalized two-body state constructed out of a state \u) 
and its time- reversed state As it is known [201 1321 [33] . the Strutinsky method averages 
the density matrix over an energy interval corresponding roughly to the distance between two 
major shells. Implicitly the same holds if the equivalent Wigner-Kirkwood expansion (TF 
approximation at lowest order) is used for obtaining p E . 

As far as we are interested in the semiclassical counterpart of the density matrix pe on the 
energy shell, we start considering its Wigner transform /g(r,p). In order to obtain the pure 
TF approximation, we differentiate with respect to E the Wigner-Kirkwood expansion of the 
full single-particle one-body density matrix p = Q(E — H) retaining only the leading term, 
which reads after normalization: 

/r(r,p) = J^(£-^), (46) 

where H c i = p 2 /2m* + V(r) is the classical Hamiltonian of independent particles with a constant 
effective mass m* moving in an external potential well. Integration over the momentum yields 
the local density on the energy shell: 

TF i \ _ 1 f j3 fTF(~ \ m*kjr;(r) 



p ^'( r ) = / d'pfViv, p) = 7,\ , (47) 
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where the local momentum at the energy E is 



v(E,E>) =/^/S/r('.P)"(P-P , )/F(',P') (50) 



Mr) = E^ = ^ (B -V(r)) (48) 

and the level density g(-E) is given by the integral of the local level density g(E, r) 

m=Jd 3 rg(E,r) = Jd\^^ (49) 

Now we proceed to calculate the average pairing matrix elements v(E,E') of the Gogny D1S 
force [M] which is known to reproduce the experimental gap values when used in microscopic 
Hartree-Fock-Bogolyubov calculations [35]. Starting from fj4"5l) and following the method ex- 
plained in detail in Ref. [3Tj . one arrives in TF approximation at 

d 3 pd 3 p' 
(2-Khy 

where f^ F is given by equation fj46|) and v (p — p') is the Fourier transform of the particle- 
particle part of the Gogny force which describes the pairing. As far as the only dependence on 
p and p' is in v (p — p'), we can average over the angle between p and p' as follows: 

v(p,p f ) = ^Jv(p-p')dQ. (51) 

This result and the fact that the TF on shell density (|4"6]) can be recast as f^, F (r, p) = m*5(p — 
Pe)/{ci{E)pe) allows one to perform easily the angular integral in ( 1501) . with the result 

1 1 /2m* \ 2 r 

v(E < E ' )= mtw^\i?) r^*^ (52) 

which in the particular case of the Gogny force reads: 

, ^ Zj 1 (2m*\ 2 fRt 2 fajkl + fc|,) . n 2 k E k E > f . 
( ' Jo drrQXp{ 4 — > ( 53 ) 

where Rt is the classical turning point and Zi = 7r 3 / 2 /^(Wi — B{ — Hi + Mi). The factors 
Zi correspond to pairing in the S = and T = 1 channel and are written in terms of the 
parameters of the Gogny force Wi, B iy Hi, M c and fa [31]. 
The semiclassical TF gap equation reads 

A(E) = - dE'v(E, E')g(E') 1 ) _ (54) 

Jv ± 2y/(E' - fi) 2 + A(E') 2 
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where V\ and V2 are the lower and upper limits of the pairing window and fi is the chemical 
potential which is obtained by the condition of the neutron (proton) number, i.e. integrating 
the corresponding TF level density up to the Fermi level: 

N T = 2 g(E)dE, (55) 

where Vo is the bottom of the potential V(R) and Ep T is the Fermi level for each type of 
particles (r =n,p). 

The pairing density in the TF approximation is given by: 

k(E) = - A(£) , (56) 
2s] (E - rt 2 + A(Ef 

Next we obtain the gap and the pairing density in coordinate space by integrating over E 
the gap and pairing density given by Eqs. ( |54l) and ( l56l) . repectively, weighted with the local 
level density g([E,R): 

A(r) = J dEg(E,r)A(E) (57) 

and 

«(r) = J dEg(E,r)K(E), (58) 

In the calculation of the energy and B (Ml) factors of the scissors mode the zeroth and 
fourth order moments of the pairing density for each kind of nucleons are needed, they read: 

k = J d 3 rk([r) =4 J dEk(E) J d 3 rg(E,r) (59) 

and 

h = J d 3 rk(r)x 2 z 2 = 4 J dER(E) J d 3 rx 2 z 2 g(E,r), (60) 

where the factor four takes into account the spin-isospin degeneracy. In order to calculate 
Eqs. f l59|) and (1601 . we use a single particle potential of harmonic oscillator type. The calculation 
of R(E) (Eq. (|56p ) is carried out in spherical symmetry and the deformation is included in 
g(E,R) in order to obtain the moments of the pairing density. With the harmonic oscillator 
potential the integral in coordinate space can be done analytically and the calculation of the 
ko and h± moments reduce to the following integrals over E: 

fco- T^A I dEE 2 K{E) (61) 
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k ± = l(^) TT^ml dEE 4 R(E), (02) 



and 

where io 2 x and uj 2 z are given by formula (1221) and the deformation S is taken from the experiment 
|17j . In all the calculation an effective mass m* = 0.8m has been used. 

Appendix D: Pair correlations in the superfluid model of 
deformed atomic nuclei 

The microscopical calculations of integrals and k\ (1191) are performed with the single particle 
wave functions of the deformed nuclei. These functions are obtained by the numerical solution 
of the Schrodinger equation with the axially deformed Woods-Saxon potential including the 
spin-orbit interaction: 

V W s{r) = -V /{1 + exp[«(r - R(6))]}, 

R(9) = R [l + /3 + (3 2 Y 20 (6) + f3 4 Y 40 (9)}, 

Vi s {r) = -K(p x a)W ws . 

Here R = r A 1//3 ; the constant /3 is introduced to ensure the volume conservation; ft 2 and (3a 
are quadrupole and hexadecapole deformation parameters. The details of the method of the 
solution can be found in [36J E7J [38l 139] . 

The wave function of the deformed nucleus is represented in a form of a superposition 

|i'>=*ft = E a 35-*Sy. ( 63 ) 

nlj 



where a n ^ are the expansion coefficients and 



= Rndr)^ (64) 

are the wave functions of the spherical basis. Here are the spherical spinors and R n ij(r) are 
eigenf unctions of the radial part of the Schrodinger equation with the spherical Woods-Saxon 
potential. This paper deals with the single particle matrix elements of the type < v\F(t)Y\q\v >, 
with F(r) and Y"ao being the radial part and the spherical function respectively. 
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Table 2: 2 , f3±, Vq (in MeV), r (in fm) and a (in fm x ) are parameters of the deformed Woods- 
Saxon potential; G (in MeV) are pairing strength constants; k are spin-orbit parameters in fm 2 . 



A 


ft 


ft 


Neutron system 


Proton system 






a 


K 


G 




^0 


a 


K 


G 


135 


0.14 





47.0 


1.26 


1.67 


0.40 


0.141 


58.0 


1.24 


1.578 


0.35 


0.150 


147 


0.13 





46.8 


1.26 


1.67 


0.40 


0.123 


57.4 


1.24 


1.578 


0.35 


0.138 


155 


0.30 


0.04 


47.2 


1.26 


1.67 


0.40 


0.115 


59.2 


1.24 


1.69 


0.36 


0.154 


165 


0.28 


0.02 


44.8 


1.26 


1.67 


0.43 


0.110 


59.2 


1.25 


1.63 


0.355 


0.132 


173 


0.26 


-0.02 


44.8 


1.26 


1.67 


0.42 


0.108 


59.2 


1.25 


1.59 


0.32 


0.133 


181 


0.20 


-0.03 


43.4 


1.26 


1.67 


0.40 


0.106 


59.8 


1.24 


1.67 


0.33 


0.130 


193 


0.14 





43.4 


1.26 


1.67 


0.40 


0.101 


59.8 


1.24 


1.67 


0.33 


0.121 



Pair correlations are taken into account in the frame of the BCS theory. The interaction 
leading to the superfluid pairing correlations acts between the particles in time-reversed con- 
jugate states. The pairing matrix element G(u+, v— ; v' — , z/+) is usually assumed to be a 
constant G independent of v and v 1 [115]. In this approximation the gap A does not depend 
on v either: A = GJ2v u u v u- Then the equations determining A and the chemical potential \i 
read: 

- = E- 

G ,, e,. 



Here N is the number of particles, u u and v v are the coefficients of the Bogoliubov canonical 
transformation, t v = (A 2 + [E v — n) 2 ) 1 ^ 2 is the quasiparticle energy, E v is the energy of the 
single particle state (63). After all transformations the integrals k$ and k\ can be written as 

kl = 4 / d 3 r<(r) = AJ2^v u J2Ki,) 2 , (67) 



nlj 



k T 4 =4j d 3 rx 2 z 2 i^(r) = 

^-Y, u » v » < H^oo + - ;koK (r)k >= 

15 ^7 7 7 

^UQpVnp E a nlj a n 'l'j'^n'l'j' E ^^jn,\0^j'l/2,\0J (^) 
Q,p nljn'l'j' A=0,2,4 
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where I*, 3 Vjl = J r 6 R n ij(r)R n >i>y(r)dr, coefficients L x are: L = ^, L 2 = ^, L 4 = -|§, 
C/™ is the Clebsh-Gordan coefficient. 

Jl " L l iJ2 lllj 2 

All considered nuclei were divided into several groups. The deformation of the potential 
Vws i n each group was chosen close to the average value of experimental deformations of nuclei 
of the group. The schemes of the single particle states for each of these groups were calculated 
with the fixed set of parameters, each set being fitted to achieve a correct sequence of the single 
particle levels in deformed nuclei. The sets of these parameters for neutron and proton systems 
are given in Table 2. The parameters for A=155, 165, 173, 181 are taken from Ref. [3D]. All 
discrete and quasidiscrete levels in the interval from the bottom of the potential well up to the 
energy 22 MeV were taken into account in the calculations. Such basis allows one to study low 
lying states and giant resonances as well. The pairing strength constants G were adjusted as a 
function of the size of the single particle basis. 
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